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We introduce a model independent method for the determination of the hadronic contribution to the 
QED running coupling, Aq:had(M|), requiring no e + e _ annihilation data as input. This is achieved 
by calculating the heavy-quark contributions entirely in perturbative QCD, whilst the light-quark 
resonance piece is determined using available lattice QCD results. Future reduction in the current 
uncertainties in the latter shall turn this method into a valuable alternative to the standard approach. 
Subsequently, we find that the precision of current determinations of AaHAD(Mf) can be improved 
by some 20% by computing the heavy-quark pieces in PQCD, whilst using e + e~ data only for the 
low-energy light-quark sector. We obtain in this case Ao:had(M§ ) = 275.7(0.8) x 10~ 4 , which 
currently is the most precise value of AajjAt>{Mz)- 

PACS numbers: 13. 40. Em, 12.20.Ds, 13. 66. Be, 13.66.Jn. 12.20.-m 



I. INTRODUCTION 



Of the subset of three parameters that enter the elec- 
troweak sector of the Standard Model (SM) of particle 
physics, Gf,Mz and a(M§), the least precisely known 
is the electromagnetic coupling at the Z boson mass, 
a(Af§). This is primarily due to hadronic contribu- 
tions which are not calculable using perturbative QCD 
(PQCD). Increasing the precision of a(M§) is impor- 
tant for, amongst other things, obtaining a Standard 
Model fit of the Higgs mass. Currently, there is a mi- 
nor tension between the recently measured mass of a po- 
tential Higgs boson, M H = 126.0(0.4)(0.4) GeV [1, or 
M H = 125.3(0.4)(0.5) GeV @, and a mass of 91+ljGeV 
(at the 68% confidence level) obtained from global SM 
fits to electroweak precision data [|[. 
The running QED coupling a can be parameterized as 

a(0) 



a(s) 



(1) 



1 - Aa L (s) - A«had(s) 

where AaL is the leptonic contribution, which can be 
determined with high precision in perturbation theory, 
and AanAD is the hadronic term. Of particular interest 
is the QED coupling at the scale Mz- Denoting a = a(0) 
in the sequel, AanAD(Af§) can be written as 

AaHAD(Mf) = Ana {n(0) - Re p(Mf)]} , (2) 

where IL(s) is the standard electromagnetic current cor- 
relator 



l J d A X e^{Q\T{j^{x),j^{Q))\Q) 

(<?m^ - q 2 g^)tt(q 2 ) , (3) 



with j^ M (x) = J^f Qff( x )7vf( x )> an d the sum is over 
all quark flavors / = {u, d, s, c, b, t}, with charges Qf. 
Invoking analyticity and unitarity for n(s), and using the 
optical theorem, i.e. R(s) = 127rImII(s), where R(s) is 
the normalized e + e~ cross-section, one can write Eq.Q 
as a dispersion integral Q 



Aa H AD(M 



a M 2 Z 
3tt 



P 



R(s) 



s{Ml - s) 



ds 



(4) 



where P denotes the principal part of the integral. This 
dispersion relation is useful as it only requires knowl- 
edge of R(s), which can be determined experimentally. 
The standard approach to determining AanAu(Mz) * s t° 
evaluate Eq.Q making use of e + e~ annihilation data for 
R(s) in the resonance regions, and either use the PQCD 
prediction for R(s) above these regions (see e.g. 0]), 
or make use of all the available e + e~ data and fill in 
the gaps using the PQCD prediction (see e.g. [E @)- 
Since the use of data is the primary source of uncer- 
tainty, other analyses have attempted to reduce the de- 
pendence of AanAD(-^i) 011 e+e ~ data by a variety of 
methods that place a greater emphasis on PQCD. One 
approach in this direction is to subtract a polynomial 
from the weight function in Eq.Q to reduce the im- 
pact of the data contribution. In order to compensate, 
this polynomial weighted integral is added to the right 
hand side of Eq.Q and evaluated in PQCD (plus non- 
perturbative corrections given in the framework of the 
Operator Product Expansion) using a circular contour 
integral (see e.g. @, HI). Another approach is to first 
calculate Acvhad(— sq) (sq > with sq large enough for 
PQCD to be valid), whose weight function deemphasizes 
the low-energy region. Subsequently, Acvhad( — so) is run 
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FIG. 1: An example of a singlet three-gluon exchange diagram 
in QCD (left), and the lowest order QED correction to the 
vacuum polarization II(s) (right). 



to Aq;had(-^§ ) using the PQCD prediction of the Adlcr 
function 0. 

The purpose of this paper is two- fold. First, to calcu- 
late the complete heavy quark (charm, bottom, and top) 
contributions to AanAD(Aff) using only PQCD, which 
to our knowledge has not been done before. Interest- 
ingly, this will significantly reduce the total uncertainty 
in Acchad{Mz ), as the use of e + e~ data in the charm- 
quark region leads to an error equivalent to that from 
the use of e + e~ data in the light-quark resonance region. 
Second, to show how existing Lattice QCD (LQCD) cal- 
culations involved in the evaluation of the hadronic con- 
tribution to g— 2 of the muon can be used to calculate the 
light-quark contribution to Aq;had(-M§) entirely from 
theory. This will allow for the first model-independent 
determination of AanAn(M^) that makes no use at all 
of e + e~ cross-section data. At present, though, current 
precision of LQCD results do not allow this method to 
compete with the standard approach. 

We begin by considering the heavy quark contribution 
to Eq.((5]), which can be written as 

Aa( I / | D (Ml)=4 7 ra(n(/)(0)-i?e[nW)(Ml)]) , (5) 

where / = {c, 6, t} are the heavy quark flavors, and in 
the sequel it should be understood that it is the real part 
of the correlators that enters in the time-like region. It 
should be noticed that one only needs knowledge of the 
correlator at s = and at s = M% . The latter scale is 
way above either the charm- or the bottom-quark pair 
production resonance region, so that one can safely use 
the high-energy expansion of the heavy quark correlator. 
This is known to O(a^) (with partial results at (D(a A s )). 
In addition, n^'(0) has also been calculated in PQCD 
to 0{a^). The other key inputs are the recent high pre- 
cision bottom- and charm-quark masses obtained from 
LQCD [Hj in the MS-scheme. This scheme will be used 
here in all PQCD calculations. Prior to these LQCD 
determinations, the charm- and bottom-quark masses 
were obtained using e + e~ data, a procedure we wish to 
avoid as we aim at an entirely theoretical determination 
of Aa>jl Ar) (Mz). As explained later, it turns out that 
Eq.([5]) is problematic for the charm-quark contribution 
(but not for the bottom- or the top-quark). The reason 
being its strong dependence on the renormalization scale, 
which must be the same for both lL c )(0) and n( c '(M|). 
Therefore, we introduce two additional approaches in the 
charm-quark sector which are significantly less sensitive 



to this problem. The first is inspired by the Adler func- 
tion approach of @, to wit. We note that 



^Aa<l D ( S ) = -4,4nW( S ) = 



where the real part is understood, and D^(s) is the 
Adler function in the charm-quark channel. Integrating 
Eq.((6]) gives 



a41 d (a/ 2 2 ) 



[A4 c | D (M z 2 )-Aa^ D ( So )] 
3tt 



a "had( s o) 



ds 



+ 4™(n (c) (o)-n (c) (s )) . (7) 

We choose s large enough so that PQCD is valid, but 
still so <C M§ . One can then use one scale for the second 
term on the right hand side above, whilst another scale 
for integrating over the Adler function (one could also use 
a running scale, e.g. /j 2 = s). The second, but similar 
approach, is to use Cauchy's residue theorem to rewrite 
the dispersion relation Eq.Q to obtain 



n( c )(o) 



Mf 



4aM§ 



ds ■ 



n< c )( s ) 



-so s(M 2 z 



s) 



3tt 



R^{s) 



s) 



■ds 



(8) 



which is only valid for sq < M§. As usual, sq will be 
taken large enough so that PQCD is valid. Once again, 
this allows for the use of more than one scale appropri- 
ate for the different regions. In addition, R( c \s) is known 
partially up to (9(oq). The final virtue of this approach 
is that it will allow a careful region- by-region comparison 
with the standard approach based on Eq.Q. The sum 
rule Eq.® is very similar to the FESR used in precision 
charm- and bottom-quark mass determinations employ- 
ing experimental data on R(s) [HI, Gil- Hence, to deter- 
mine Aq;^2d(-^z ) entirely from theory it is essential to 
use a non-QCD sum rule determination of the charm- and 
bottom-quark masses, such as e.g. that from LQCD. The 
procedure just outlined for determining Aaj( AD (M 2 ) is 
not necessary for the bottom- and top-quark counter- 
parts, as Eq.© for / = b,t gives results that are essen- 
tially renormalization scale independent. However, we 
have checked that the FESR and the Adler function ap- 
proaches give the same result as that using Eq.© for 
the charm-quark contribution, although the latter has a 
much larger error. 



II. VECTOR CURRENT CORRELATOR IN 
QCD 



We provide in this section a summary of the available 
theoretical information on the vector current correlator 
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where 




10 11 



FIG. 2: The difference between the order O(al) and the order 
(D(a 2 s ) results for 47ra(II (c) (0) - Re[II (c) (Ml)]), i.e. the trun- 
cation error Atr, as a function of the renormalization scale 
fi. 



n (n) W = E 



n 



(») 



(12) 
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in QCD. The flavor /-quark current correlator can be 
split as 

n (/) (s) = ii^cd W + n$(*) + n<g D ( s ) , (9) 

where / G {uds, c, b, i\, ^pqcd is the PQCD contri- 
bution, II^p is the contribution from non-perturbative 
power corrections given in the framework of the Opera- 
tor Product Expansion (OPE), and n^, D are QED cor- 
rections as shown in Fig. [T] The dominant contribution 
to Eq.© is the perturbative part np^ CD . In the mass- 
less case, appropriate for the u, d, s quarks, the high en- 
ergy limit of n p ^ CD is known exactly up to order O(ot 3 s ), 
and up to a real constant to order 0(a 4 ) (th e full re- 
sult is given in computer readable form in [13|). In the 
heavy-quark case, though, one needs both the low- and 
the high-energy expansions of the correlator. In the low- 
energy limit, with a single heavy quark and nf active 
flavors, the vector correlator can be written as 



n/00 



3Q 2 



16 7T 



i=0 



Ci 



(10) 



where z = s/(4m 2 ), and to/ is the mass of the quark of 
flavor-/ in the MS scheme at the scale /i. The coefficients 
Co and C\ were determined up to 0{aP s ) in 0, 13: ^2 
in [HI , and C3 in pj} . In the high-energy limit the heavy 
quark correlator is written as the massless one with added 
quark-mass corrections 



n( s ) = q 



nW( f 



(11) 



FIG. 3: The value of n (uds) (0) in LQCD in the MS-bar 
scheme at [i = 2 GeV for different values of the pseudoscalar 
mass. Solid line is a linear fit extrapolation to the physical 
pion mass (indicated by the vertical dashed line). The errors 
from the LQCD fit parameters are added in quadrature. 



The terms n(°) and IlW above are known exactly, whilst 
n^ 2 ' is for all practical purposes also known exactly, i.e. 
the mass corrections up to 0(m 60 ) are given in (Til, [l8l - 
[20( 1 . At order 0(al), there is only partial information 

(3) (3) 

on the mass corrections. The functions Tl Q and H\ ' 
are known exactly [2l[, whilst the logarithmic terms of 



n^ 3 -* are given in (22[. The constant term has been es- 
timated using Pade approximants (23[, but it was found 
to be negligible. At order 0(a*) the logarithmic terms 



1 n 4) and n 1 



are known [2J, |25| , but not the constant 
terms. Hence, for the heavy quarks no 0{a A s ) terms will 



T (4) 



and 



be included in any contour integral. However, H 

n^ 4 ' will be included when integrating over Im(n(s)), as 
the unknown constant terms do not contribute in this 
case. At 0(a? s ), there are singlet diagrams contribut- 
ing for the first time (Fig. [T]). The logarithmic terms 
are known in the high energy case (22| . and we will use 
these in Eq.(|3]) to estimate the singlet contribution. An 

example of a lowest order contribution to IlQg D (s) is 
shown in Fig. [TJ after substituting gluons by photons. 
Finally, the leading-order non-perturbative contribution 
to n^jp(s) for heavy quarks is from the gluon condensate 
((a s /7r)G 2 ). This has been determined from data on r- 
decays [26j , and with a very large uncertainty from data 
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Contributions to Aqhad (in 


units of 10~ 4 ) 








[U O. / Vj»C V Juds 


[O.t U.O\jii\\ u dsc [y.O <±U V_re V J udscb 


f^> 4f) Hp VI j i 

L^> ^iU Vj.c V J udscb 


1 ^ 4:U Vjrc V J t 


Total 


Standard approach [3] 


79.29 ± 0.69 


60.21 ±0.51 93.50 ±0.16 


42.70 ±0.06 


-0.72 ±0.01 


275.0 ±1.0 


This work 


79.39 ± 0.68 


60.46 ±0.33 93.82 ±0.14 


42.76 ± 0.06 


-0.76 ±0.03 


275.7 ±0.8 



TABLE I: The contributions to Aqhad from different regions using either Eq.Q (standard approach), or Eq.(|8} (this work) 
which requires data only for s < 1.8 GeV 2 and PQCD above this energy. Our results in columns three and four are obtained 
using Eq.Q also for the bottom-quark contribution, with t/sq = 40 GeV. The total error takes into account the correlations 
of the uncertainties from the different regions. 



on e + e annihilation into hadrons [27] . The conservative 
value <(a s /7r)G 2 ) = (0.006 ± 0.012) GeV 4 will be used in 
the sequel. To estimate the error arising from the incom- 
plete knowledge of the correlator in PQCD (truncation 
error), we take the difference between the C( a s) an d 
C(a" _1 ) results, where n is the highest available order. 
We will check this by also varying the scale /j,. Finally, 
many errors will be 100% correlated or anti-correlated 
between different regions, such as e.g. the error in a s . 
As input, the PDG value of the Z-mass will be used, i.e. 
Mz = 91.1876(21) [28]. For the strong coupling we use 
the result from the determination of Davier et al. 0] 



a a (M|) = 0.1193(28) , 



(13) 



in order to facilitate the comparison of the final result 
for Ao!had(-^I)- This value has a far more conservative 
error than the PDG value of a s {M 2 z ) = 0.1184(7) [Hj]. 
However, this will hardly matter as our results turn out 
to be largely independent of the strong coupling. For the 
charm- and bottom-quark masses, we use the most recent 
LQCD determination Q3, m, 4) (3GeV) = 0.986(6) GeV 

and tt4 5) (10 GeV) = 3.617(25) GeV. These values are 
in very good agreement with QCD sum rule deter- 
minations [ll|, 0, Ull- For the top-quark mass, we 
use m[ 5 \m t ) = 160.0(3.5) GeV (3(j. If one employs 
either the Adler function or the FESR approach, the 
coupling and quark masses need to be run across flavor 
thresholds. This will be done using the Mathematica 
program RunDec [3lj . 



III. CHARM-QUARK CONTRIBUTION 



We consider first the evaluation of Aa^ D (Mf ) directly 
from Eq.([2]) . As mentioned previously, there is a prob- 
lem with the strong dependence of Aa^ D (M|) on the 
renormalization scale \x. In fact, if one varies [i in the 
interval [i = 2 GeV - M z then Aaj£' AT) (M z ) changes by 
0.77 x 10 -4 , which is a large variation in the context 
of the current precision of Ao;had(-^I)- We choose in- 
stead a value of fi leading to a good convergence of the 



perturbative series. Figure [21 shows the difference be- 
tween the order 0(al) and the order O(oj.) results for 
47ra(n( c '(0) - Re[n>)(Af§)]), i.e. the truncation error 
Atr, as a function of the renormalization scale (i. A value 
of fj. ~ 3 GeV ensures good convergence. Using rif = 5, 
this leads to 

a4ad(^I) = 4mfn( £ »(o)- [nW(Mi)] 



= (79.19 ±0.13 Atr ±0.03 Ac[a 

± 0.01 A ( G 2> ± O.llAmJ x 10~ 4 , (14) 

where the errors are due to truncation (Atr), and to un- 
certainties in a s (Aa s ), in the gluon condensate (A(G 2 )), 
and in the charm-quark mass (Am c ). Interestingly, this 
value corresponds to a global minimum of Aa^ AB (M z ) 
as a function of /i. In any case, given this strong 
/i-dependence we shall not use this method to deter- 
mine the charm-quark contribution. We consider in- 
stead the Adler function approach using Eq.® which 
involves two terms, a high energy and a low energy con- 
tribution. Starting with the low energy part, and using 
rif = 4, s = (9.3 GeV) 2 (which is below the bottom- 
quark threshold) , and an initial value fi = 5 GeV (to be 
made to vary in a wide range later) , we find 

47ra(n (c) (0) - n (c) (s )) = (29.57 ± 0.25 A tr 
±0.05 AQs ± 0.01 A(G 2) ± 0.12 Al - n J x 10~ 4 . (15) 

Allowing now [i to vary in the wide range fj, = (2 — 
9.3) GeV, results in a variation of the central value above 
in the range [29.31 — 29.62] x 10 -4 , which is within the 
truncation error. For the high energy term, we use rif = 5 
and an initial value /i = \M Z to obtain 



o 

3^ 



-ds = (49.91 ± 0.05 Atr ± 0.04 Aaa 
T 0.01 Afl J x 10~ 4 . (16) 



Varying fi in the interval fi = 9.3 GeV — Mz, produces a 
variation of the central value above in the range [49.82 — 
49.95] x 10~ 4 . Adding both contributions, Eqs.([T5]) and 
(fTTJl) , the total result from this approach is 



Aa 



had(^I) = (79.49 ±0.30 Atr ±0.09 Ac[a 

± 0.01 A(G 2) ±0.1lAmJ x 10~ 4 ,(17) 
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where we notice that the errors Am c in Eqs. (|15[) and 
(fT6")) are anti-correlated. The third method consists in 
using the FESR Eq.flEJl. The values n/ = 4, and /i = 
5 GeV, have been used both in the contour integral, with 
s = (9.3 GeV) 2 , as well as in the residue 11^(0). For the 
integral involving R( c \s) we use nj = 5 (the contribution 
above the top threshold is not numerically important), 
and obtain 

a41d(^I) = (79-34 ± 0.26 Atr ± 0.04 AQs 

±0.01 A<G 2) ± 0.11a™.) x 10" 4 .(18) 

Varying fi in the respective regions as done previously 
produces changes in the central value of Aa^ D (M§) 
well within the truncation error. This result will be 
adopted for the charm-quark contribution as it is has 
a much smaller uncertainty due to a s , and a slightly 
smaller truncation error. It should be mentioned that 
at s — the low energy expansion, Eq.(JTU|), is a well 
convergent power series expansion in the strong coupling. 
In the vicinity of s = 0, Eq.(fT0|) also converges well, even 
if the charm-quark is at the borderline between light 
and heavy quarks. This is a standard procedure in the 
determinations of the charm-quark mass from QCD sum 
rules or from LQCD. 



IV. BOTTOM- AND TOP-QUARK 
CONTRIBUTIONS 



In this case all three methods give essentially the same 
result for AajJ AD (M§). Using e.g. Eq.([2|) with n/ = 5 
and \x = 10 GeV gives 

a4 6 1d(^I) = 4na(nW(0)- [nW(Mf)]) 
= (12.79 ± 0.06 Atr ± 0.009 Acts 
± 0.03 Afft J x 10- 4 . (19) 

Varying fi in the very wide range \x = 10 GeV — 10Mz 
only changes this result by 0.04 x 10 -4 , which shows a 
remarkable scale independence. 

For the top quark, one can use the low energy expansion 
to calculate both 11^(0) and IlW(M|). Up to 0(a s ) the 
full analytic correlator is known. We have verified that 
there is no appreciable difference between results using 
the low-energy expansion of the correlator or using the 
full expression up to this order to determine IT^Mf ). 
At higher orders, one can reconstruct the full analytic 
behavior of the correlator using Pade approximants, as 
done in [23| at order O(a^). Using these results we find 
that it is perfectly safe to use the low energy expansion 
of the correlator. With fi = rh t and n/ = 6, we find 

AagL(Mf) = 47ra(nW(0)-IlW(M!)) 

= (-0.76 ± 0.03 Am J x 10~ 4 , (20) 



where only the uncertainty in the top-quark mass pro- 
duces a non- negligible uncertainty in Aa{j f ^ D (Af|). 



V. LIGHT-QUARK CONTRIBUTION 



In contrast to the heavy quark contributions, one can- 
not use PQCD to determine the light-quark correlator at 
low energies. There are two approaches to achieve this, 
i.e. using e + e~ data for R(s), or LQCD determinations 
of II(s) (in the space- like region), with both being used 
in the sequel. For the e + e~ data approach, and below 
the onset of PQCD, we use the integrated result of Q 
to avoid the complicated task of dealing with the vast 
amount of e + e~ data available (for a recent independent 
analysis see [13 )■ The result of @, integrated up to the 
PQCD threshold y/s = 1.8 GeV is 

«MS / (1 - 8GcV)2 fl&ffOO d: 

Sir J s(Mf-s) 

= (55.02 ±0.66) x 10" 4 . (21) 

Above the PQCD threshold, and using the massless order 
C(a 4 ) PQCD expression for R^ uds ^{s) we find 

aM% /"» 4qcd(*) ^ 

37r Al.SGoV) 2 S ( M Z ~ S ) 

= (129.26 ± 0.16 Atr ± 0.29 Aq J x 10 -4 , (22) 

which added to Eq. (|2T|) gives the total light-quark con- 
tribution 

a «had( m !) = (184.28 ±0.66 data ±0.16 Atr 

± 0.29 Aa.) x 10" 4 . (23) 

A large effort is currently underway to determine 
n(«<* s )( s ) in the space-like region using LQCD. A key 
aim is to provide a first-principles determination of 
the hadronic contribution to the g — 2 of the muon 
[3214361 ]. Here we describe two methods for obtaining 
Aa^Q (-M|), entirely from theory, from a combination 
of LQCD and PQCD. This is inspired by [38|, where an 
entirely theoretical determination of g — 2 was proposed. 
It must be emphasized, though, that LQCD results are 
currently not precise enough to compete with the e + e~ 
approach. For instance, one source of uncertainty arises 
from disconnected Feynman diagrams, which are cur- 
rently not included in LQCD calculations, and which lead 
to an estimated 10% systematic uncertainty (34|. 
The first method is based on the FESR Eq.@, with 
n(" ds )(0) determined from LQCD, and the two integrals 
computed in PQCD. In LQCD it is not possible to cal- 
culate directly I\^ uds ^ (0) . Instead, the correlator is com- 
puted for values very close to s = 0, and then these 
results are fitted and extrapolated to the origin to obtain 



G 



Authors 


(M|) (in units of 10~ 4 ) Method 


Groote et al. (199o) [YJ 


277.6(4.1) 


r\ /—\ /-^-r\ ti / ■ "I • l ■ 1 j 11111 l*i \ 

PQCD driven (using polynomial weights and global duality) 


Kuhn et al. (1998) \oj 


277.5(1.7) 


r»/^/nFv i" / ■ l "i "ii iiiii i*j \ 

PQCD driven (using polynomial weights and global duality) 


Burkhardt et al. (2U11) 39 


275.0(3.3) 


Data driven 


Irocomz et ai. (2005) |40| 


274.9(1.2) 


Data driven 


Jegerlehner (2008) [9J 


275.94(2.19) 


Data driven 


Jegerlehner (2011) 1411 


274 QSfl Sfil 


A H Ipt Ti i n i nn 3nnrna pn 

iiUlcl ILIIILjLUJII CIUU1UCIL.11 


Hagiwara et al. (2011) [5] 


276.26(1.38) 


Data driven 


Davier et al. (2011) [3] 


275.7(1.0) 


PQCD m range 1.8 < ys < 3.7Gev, and for > 5 GeV, 






otherwise data. 


This work 


276.5(0.8) 


Data for < 1.8 GeV, and PQCD for \/s > 1.8 GeV. 


This work 


~ 273 


LQCD + PQCD 


TABLE II: Some five-flavor 


results of previous analyses 


of Ao?had(M§) together with our determinations. Only the latest 



available results from each collaboration are quoted. No error is given in the last line (LQCD + PQCD) in view of the large 
uncertainties from LQCD. 



n u d s (0). For instance, the phenomenologically inspired 
fitting function used in j34j is of the form 



with fit parameters given in [341 ] . An alternative, model 
independent approach to extrapolating LQCD data is 
based on Pade approximants [37| . This approach would 
be appropriate in future precision determinations based 
on improved LQCD data. An important observation is 
that neither H u ds(0) nor the contour integral in Eq.© 
are observable quantities. Therefore, it is essential to 
compute both of these quantities in the same rcnormal- 
ization scheme, and at the same scale, so that the observ- 
able difference between H u d s (0) and the contour integral 
is scheme-independent. A problem arises because PQCD 
schemes, such as MS, are not easy to relate to LQCD 
rcnormalization schemes. The latter lead to a predic- 
tion of U( uds >(s) which differs from the MS results by 
the constant LV" ds )(0), which is precisely what is needed 
in Eq.(j8]). The standard approach to fix this constant 
is to impose agreement between PQCD and LQCD re- 
sults at some value s = — s* where PQCD is expected to 
be valid. This procedure would then allow for a deter- 
mination of ]T> ds )(0) from LQCD. Above s ~ -2 GeV 2 
LQCD results already are in agreement with PQCD, so 
to determine ±l( utis )(0) we choose s* = -3.5 GeV 2 to be 



on the safe side, together with s = (3.72 GeV) 2 which 
corresponds to the onset of the charm-quark region, and 
n f = 3. The renormalization scale was varied in the wide 
range between the r-lepton mass and the charm thresh- 
old, i.e. yu = 1.77— 3.7 GeV. This produces a negligible 
change in Aq^' (Mi) of 0.05 x 10" 4 . There are eight 
available LQCD fits [3J] for different pion masses, with 
each giving a result for IT" ds )(0) in need of extrapolation 
to the actual physical value of the pion mass. Using a 
simple linear extrapolation fit gives the results shown in 
Fig. [3l leading to 

n< uds )(0) ~ 0.08758 Gu = 2 GeV) . (25) 

Using this value in Eq.©, together with a PQCD evalu- 
ation of the integrals, gives 

Aa^(M|) ~ 181 x 10~ 4 . (26) 

No error is given above in view of the current uncertain- 
ties in LQCD. 

The second method to determine Aa^j (Mf) entirely 
from theory makes use of the Adler function. One im- 
portant advantage over the FESR approach is that there 
is no need to enforce the matching between PQCD and 
LQCD at any point. In order to involve the Adler func- 
tion we write Aa^^(M^) as 



Aa^SCMf) A<d } (-*,)+ [A 



(uds 



Ana 



n 



(uds) 
LQCD 



(uds) 
*HAD 



(o)-niS(-«o) 



/ \ a (uds). 



LQCD \ 



3- 



■-o)] + W 

Juds 



D 



PQCD 



(a) 



ds 



{Ml 



3tt 



Act 



Ml 



(uds), v1 
HAD y S 0)\ 
(uds 



D 



pqcd( s ) 



ds, 



(27) 
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where Sq is large enough for PQCD be valid, and the real 
part of the expression in the last line is to be understood. 
Notice that the line integral in the interval (— sq, +so) is 
well defined. In fact, since D(s) is an analytic function 
the integration was performed on a semi-circular contour 
of radius |so|, avoiding the origin. Evaluating each of 
the three terms on the right hand side of Eq. ([2~T|) . with 



the gluon condensate, are negligible. This result can 
be compared with Aoihad(M§) = (275.0 ± 1.0) x 10~ 4 
from [3| (for nt = 6) (other five quark flavor results are 
listed in Table II). The primary reason for this 20% 
reduction in uncertainty is our PQCD calculation of the 
contribution of the charm-quark resonance region, which 
is given in Table |U 



so 



= -3.5 GeV, we find 

natoi-nS'f-.o) 



~ 52.32 x 10 -4 (28) We comment in closing on the relation between 

a I a »9 \ i,i i r .i ti- mi. 



<0 D 



(uds) 



(s) 



a_ I -PQC D V^ & = 2 56 xl0 - 4) (2g) 



a 
3tt 



ds = 125.95 x 10" 4 , (30) 



which add up to 

A<J 



(31) 



as already given in Eq. (|26|1 . The errors in the total 
PQCD contributions are 



Q 

3^ 



n (uds) , . 
f^PQCDlfJ ds Q_ 

S 37T 



(128.5 ± 0.2 Atr ± 0.3a«.) x 10" 



ds 



(32) 



The contribution of the gluon condensate is at the level 
of one order of magnitude smaller than the uncertainty 
in a s , hence it can be neglected. Once LQCD deter- 
minations of Eq. (|28[) achieve enough accuracy, it would 
become possible to determine AanAD(M z ) entirely from 
theory, after adding to the LQCD light-quark contribu- 
tion the heavy-quark results Eqs. ([r5]l . COB , and (poll. 



Ao;had(-M§) and the value of the Higgs mass. To be- 
gin with, as is well known both AaHAD(-^jf§) ano - the 
QCD strong coupling a s (M z ) enter into the global SM 
fit to electroweak precision data. At first sight one could 
suspect that the central role played by PQCD in our 
approach could lead to a stronger correlation between 
Ao;had(-M§) and a s (M z ). However, this is not the case. 
In fact, quite the contrary, i.e. in comparison with the 
standard approach our FESR method reduces the over- 
all dependence of Acchad(M z ) on a s (M z ). This hap- 
pens because the a s (Mz) dependence of the combina- 
tion of contour integral plus Cauchy residue in Eq.® 
is anti-correlated with the a s (Mz) dependence of the 
integral involving R(s). Hence, there is some cancel- 
lation of the a s (Mz) dependence, which does not take 
place in the standard approach which uses data to de- 
termine the heavy-quark resonance contribution. Quan- 
titatively, the functional dependence of the central value 
of Ao;had(jW§) on the value of a s (Mz) in the standard 
approach [3| is approximately 132.1[a s (Mz) — 0.1193] x 
10~ 4 , while in our method it becomes approximately 
107.7[a s (M z ) - 0.1193] x 10" 4 . 

The correlation between Aa^ D (il/§) and the logarithm 
of the Higgs mass, In My, was found in [42j to be —0.395 
(for an earlier determination see [HI). Using our result 

for Aaj?l D (M|) in Table II would lead to a Higgs mass 
Mjj ~ 87 GeV, somewhat lower than the value from Q 
Mh — 91^23 GeV, thus increasing the tension between 
a possible Higgs of mas Mh — 126 GeV, and the fitted 
Higgs mass. 



VI. CONCLUSIONS 



Adding up all of the contributions, i.e. Eas. (fT8l) . (|T9|) . 
pO)) . and d23j , the final result for Aa HAD (M z ) is 

Aa H AD(Mf) = (275.7 ±0.66 A data±0.44 Atr ±0.26 AQ 
± 0.11 A „J x 10~ 4 

= (275.7 ±0.8) x 10~ 4 , (33) 

where rif = 6 has been used, and the uncertainties due 
to the bottom- and the top-quark masses, and due to 
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